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Abstract 



r-| \ The two-dimensional M = (2, 2) Wess-Zumino (WZ) model with a cubic 

superpotential is numerically studied with a momentum-cutoff regularization 
'nI" ■ that preserves supersymmetry. A numerical algorithm based on the Nicolai 

map is employed and the resulting configurations have no autocorrelation. 



'\0 ', This system is believed to flow to an A/" = (2, 2) superconformal field theory 

(SCFT) in the infrared (IR), the A2 model. From a finite-size scaling analysis 
of the susceptibility of the scalar field in the WZ model, we determine 1 — 

Q ■ h — h = 0.616(25) (13) for the conformal dimensions h and /i, while 1 — 

h — h = 0.666 . . . for the A2 model. We also measure the central charge in 
the IR region from a correlation function between conserved supercurrents 
and obtain c = 1.09(14)(31) (c = 1 for the A2 model). These results are 

rS ' consistent with the conjectured emergence of the A2 model, and at the same 

c^ ■ time demonstrate that numerical studies can be complementary to analytical 

investigations for this two-dimensional supersymmetric field theory. 
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1. Introduction 



It is believed that the infrared (IR) hmit of the two-dimensional M = 
(2, 2) Wess-Zumino mode0 (2D J\f = (2, 2) WZ model) with a quasi-ho- 
mogeneous superpotentia iis a non-trivial Af = (2,2) superconformal field 
theory (SOFT) 0, S, i, STB, 0, B H [ij. See Section 19.4 of Ref. [nj and 



Section 14.4 of Ref. [l2| for reviews. This Landau-Ginzburg (LG) descrip- 
tion [l3[ of A/" = (2, 2) SCFT is a remarkable non-perturbative phenomenon 
in field theory and physically, for example, provides a basis for application 
of the gauged linear sigma model [14] to the Calabi-Yau compactification. 
Although the emergence of SCFT has been tested in various ways, it is 
very difficult to confirm this phenomenon directly in correlation functions, 
because the 2D WZ model is strongly coupled in low energies. Applica- 
tion of conventional numerical techniques (such as the lattice) would not be 
straightforward either, because supersymmetry (SUSY) must be essential in 
the above non-perturbative dynamics. 

In a recent interesting paper [15|], Kawai and Kikukawa revisited this 
problem and they computed non-perturbatively some correlation functions 



in the 2D WZ model by employing a lattice formulation of Ref. |16[. They 



considered the 2D A/" = (2, 2) WZ model with a massless cubic superpotential 

]¥{<!>) = ^$3, (1.1) 

which, according to the conjectured correspondence, should provide a LG 
description of a pair of the A^ = 2 c = 1 minimal models, where one is 
left- moving and the other is right-moving (the so-called A2 model). In the 
IR limit, the scalar field in the WZ model A is identified with a chiral pri- 
mary field in the A2 model with the conformal dimensions {h, h) = (1/6, 1/6) 
and U{1) charges (g, g) = (1/3,1/3). (The complex conjugate A* is identi- 
fied with an anti-chiral primary field with {h, h) = (1/6, 1/6) and (g, q) = 
(—1/3, —1/3).) The authors of Ref. [l5| obtained finite-size scalings of scalar 
two-point functions which are remarkably consistent with the above SCFT 
correspondence, thus demonstrated the power of a lattice formulation of this 



^This system is obtained by a dimensional reduction of the four-dimensional Wess- 
Zumino model [Ij] from four dimensions to two dimensions. 

^A polynomial W^(<&) of variables ^i {I — 1, 2, ... , N) is called quasi-homogeneous 
when there exist some weights cuj such that W{^i — > A"^<i>/) = AW{^). 



supersymmetric field theorylfl 

In this paper, motivated by the success of Ref. Ig], we study the 2D 
M = (2, 2) WZ model with massless cubic superpotential (11.11) numerically. 
We employ a non-perturbative formulation advocated in Ref. [251] that uses a 
simple momentum cutoff regularization. Although there is an issue concern- 
ing the locality in this formulation, the restoration of an ex pec ted locality 
property can be shown at least within perturbation theory [25|. This for- 
mulation possesses very nice symmetry properties: it exactly preserves full 
SUSY, translational invariance, and linear internal symmetries such as the 
i?-symmetry. We believe that these nice symmetry properties are especially 
useful in defining Noether currents in the regularized framework. In fact, 
by defining conserved supercurrents and identifying the component of the 
superconformal currents in the IR limit, we numerically measure the central 
charge of the system in the IR region: Together with a measurement of the 
conformal dimension, this forms a main result of the present paper. 

Throughout this paper, Greek indices from the middle of the alphabet, 
fi, u, . . . run over and 1. Greek indices from the beginning a, /3, ... are 
for spinor indices and run over 1 and 2. Repeated indices are not summed 
over unless explicit summation symbol is indicated. We extensively use the 
complex coordinates defined by 

z = xq + ixi, z = xo — ixi, (1.2) 

and 

d, = ^{do-td,), d, = ^ido + td,). (1.3) 

Conjugate momenta are defined by 

Pz = ^{po-ipi), Pz = -{po + ipi)- (1-4) 

Two-dimensional gamma matrices are defined by 

70 =r J), 71 =r, o), (1.5) 



■^For preceding numerical simulations of the 2D J\f = (2, 2) WZ model with a massive 
cubic superpotential VF($) = m^'^/2 + X<^^/3, see Refs. [iS, [ll, [3 0, [El, [53] . See also 



Refs. [23j, |2J] for theoretical background. 



and 

7.= (o J), 7. -(J o)- (1.6) 

2. Supersymmetric formulation of the 2D A/" = (2, 2) WZ model 

We start by recapitulating the formulation of Ref . [Zg] . We suppose that 
the system is defined in a finite box with a physical size LqX Li. The Fourier 
modes f{p) of a periodic function in the box are defined by 

^(^) = 7VE^'''/(p)' /(p) = fd'xe-'P^fix), (2.1) 

where the momentum p takes discrete values 

Pm = 7^^m, n^ = 0,±l,±2,..., (2.2) 

and by the definition, 

f*ip) = h-p)\ (2.3) 

where the left-hand side denotes the Fourier transformation of the complex 
conjugate of f{x), f{x)*. In the present formulation ^25|, we restrict the 
momentum p by an ultraviolet (UV) cutoff A, 

- A < Pm < A, for /i = 0, 1. (2.4) 

We parametrize this UV cutoff by a "lattice spacing" a, 

A ^ -, (2.5) 

a 

although we do not assume an underlying spacetime lattice structure in this 
paper (see below). Throughout this paper, all dimensionful quantities are 
measured in units of the lattice spacing. In particular, if we define the number 
of lattice points A^^ by 

L^ = N^a, (2.6) 

then imposing Eq. (12. 4p implies 

-^<n,<^, for /x = 0, 1, (2.7) 



and thus the number of points in the momentum grid is given by ^ 1 = 
(iVo + l)(iVi + l). 

In the present formulation of the 2D M = (2, 2) WZ model, the partition 
function is defined by 



Z 



n 

-A<p^<A 



dA{p) dA*{p) Yl d^a{p) n d^o^iP) dF{p) dF*{p) 



a=l 



o=l 



(2.8) 

where A, {ipa, "ipa), and F are the scalar, fermion and auxiliary fields, respec- 
tively. In this expression, the action S is simply the action of the continuum 
WZ model in terms of Fourier modes: 



S 



^E 



LqLi 



4:PzA*{-p)p^A{p) - F*{~p)F{p) 
- F*{-p) * W'{Ay{p) - F{-p) * W'{A){p) 



, (2.9) 



where the holomorphic function iy(y4) is the superpotential and * denotes 
the convolution 



{if I * If 2) (p) = -j—r Yl '^i(^)^2(p - q)- 



LqLi 



(2.10) 



Products contained in ly' (A) and W"{A) are also understood as convolutions. 

Since action 02. 9p is identical to that in the continuum theory, the reg- 
ularized theory (12. 8 p is manifestly invariant under all symmetries that are 
consistent with the momentum restriction (12. 4p . This is the case for symme- 
try transformations that act linearly on field variables. Thus, in the present 
formulation, SUSY, translational invariance, and the /^-symmetry (if it ex- 
ists) are exactly preserved. One can derive Ward-Takahashi (WT) identities 
associated with these symmetries in a regularized framework. 

What is sacrificed in the present formulation, on the other hand, is lo- 
cality. One sees that the kinetic terms and the interaction terms are quite 
non-local in the configuration space when the UV cutoff A is finite. In fact, 
when both integers A^o aiid Ni are odd, the present formulation is nothing 
but a two-dimensional version of a lattice formulation of the four-dimensional 



WZ model studied in Ref. 26 that is based on the SLAC derivative 12 



A detailed analysis of a one- dimensional version 
can be found in Ref. 



29 



i.e., quantum mechanics) 
Thus there is an issue of whether the present 
formulation reproduces an IR physics expected in the original target theory 
(i.e., whether it belongs to the same universality class or not). Although 
one can show 25|] that within perturbation theory the expected locality is 
restored in the limit A — t- oo for two- and three-dimensional WZ models, 
the validity of the present formulation at the non-perturbative level is not 
obvious a priori. We believe that our results in this paper will provide an 
affirmative indication regarding this question of locality. 

As a side remark, we note that the present formulation cannot be regarded 
as a spacetime lattice formulation when either Nq or Ni is even. Fourier 
modes f{p) of a function f{x) defined on a spacetime lattice is periodic in 
the Brillouin zone f{p + (7i/a)fi) = f{p), where fi denotes a unit vector in 
the n direction. However, the combination p^fXp) appearing in the action 
breaks this periodicity at the boundary of the Brillouin zone and this cannot 
be regarded as a Fourier transformation of a function defined on a lattice. 
(When Nfj^ are odd, there is no Fourier mode on the boundary of the Brillouin 
zone.) Throughout this paper we set iV^ even and this means that we lose 
a connection with a spacetime lattice formulation; we have to regard our 
configuration space as continuous (with a limited resolution). Still, Eq. 02. 8p 
provides a regularized partition function and can serve as a starting point 
for non-perturbative study. 



3. Simulation algorithm based on the Nicolai map 

After integrating over the auxiliary fields F and F* in Eq. (12. 8p . the 
partition function becomes 



Z 



n 

-A<p^<A 



dA{p) dA*{p) JJ #„(p) JJ ci4(p) 



where 



^"Z^S 



N*{-p)N{p) 



a=l 



q:=1 



(3.1) 



.(....)(-.)i^-,-lf) Ik.) 



(3.2) 



In this expression, N{p) is a function that specifies the Nicolai map 30l. l3ll. 



Nip)=2tp,A{p) + W'{Ay{p), 

N*{-p) = N{py = -2ip,A*{-p) + W\A){-p). 



(3.3) 
(3.4) 



For example, for cubic superpotential (11.11) . the exphcit form of the func- 
tion N{p) is given by 



Nip) = iipo - ipi)A{p) + A—— J2 A*iq)A*{p - q) 



LqLi 



{po - tpi)A{p) + X—— J2 MqTM-p - qT- 



LoL 



(3.5) 



We then note 
1 



S 



LqLi 



E 



dNjp) dNjp) ,^, 

N*{-p)Nip) + {^,M{-P) I ^ 1% 2 1 (^) 

dA(p) dA*{p) / \^2 



and therefore, after integrating over fermion fields 

1 



Z = 



Yl \dA{p)dA*{p) 



-A<p^,<A 



exp 



LqLi 



J2N*i-p)Nip) 



det 



(3.6) 

d{N,N*) 



d{A,A*) 

(3-7) 
In this partition function, we may change integration variables from {A{p), A*{p)) 

to {N{p), N*{p)). Then the Jacobian associated with this change of variables 

precisely cancels the absolute value of the fermion determinant. In this way, 

we arrive at the expression 



Z 



Yl \dN{p)dN*{p) 

-A<Pp<A 



exp 



LqLi 



J^NipTNip) 



X 



2_] sign det 



d{N,N*) 



d{A,A*) 

where A{p)i and A*{p)i (z = 1, 2, . . . ) are solutions of 

2ip,A{p) + W'{A){-py - N{p) = 
-2ip,A{py + W'{A){-p) - N{py = 



, (3.8) 



jT. -iT-i t-^ -"T- i 



(3.9) 
(3.10) 



Representation fl3.8p for the partition function (and a similar representa- 
tion for expectation values of observables) gives rise to the following simula- 
tion algorithm 171. Il5|: 



(i) Generate Gaussian random numbers {N{p),N{p)*). 



(ii) Find (numerically) all the solutions {A{p)i, A{p)*) (i = 1, 2, . . . ) of Eqs. (13. 9p 
and dXTU]) . 



(iii) Then compute the sums 



y signdet 



d{N,N*) 



and 



y signdet 



d{A,A 

d{N,N* 



A=Ai,A*=A* 



d{A,A* 
where 0{A, A*) is an observable. 



0{A,A* 



(3.11) 



(3.12) 



A=Ai,A*=A* 



(iv) Repeat the steps from (i) and take an average over the Gaussian random 
numbers. The expectation value of an observable O is given by 



(O) 



,,E. signdet 'i^OiA A*] 



J\ — jt.^,jt. — J\ i 



N 



(3.13) 



where {■)n denotes the average over the Gaussian random numbers and 



A 



Y,i sign det 



d{N,N*) 



d{A,A*) 



J\. — J±l^J\. — J\ i 



N 



(1) 



(3.14) 



N 



As Eq. (I3.13P shows, in the present algorithm the expectation value of an ob- 
servable (O) is given by the ratio of two expectation values with respect to 
the Gaussian random numbers. In this paper, we apply a simple propagation 
of error rule to the ratio to estimate the statistical error in (O). This proce- 
dure must be good enough because the statistical error in the denominator 
of the ratio. A, is generally quite small (see below). 



8 



Since the above simulation algorithm is based on the generation of the 
Gaussian random numbers, there is no autocorrelation among generated con- 
figurations and there is no critical slowing down; this is an overwhelming ad- 
vantage of the present simulation algorithm^ Another interesting feature of 
this algorithm is that the "normalized partition function" A in Eq. fl3.14p can 
be argued to give the Witten index tr(— 1)-'^ |35|o For the massive free the- 
ory iy($) = 'm$^/2, one easily sees that A = 1 and reproduces the correct 
Witten index tr(— l)''^ = 1. Thus assuming that the proportionality constant 
between the partition function defined by the functional integral and tr(— 1)^ 
is independent of the interaction, we have A = tr(— 1)^. It is known that for 
a homogeneous superpotential iy($) = A$"/n, tr(— 1)'^ = n — 1. Therefore, 



the relation A = n — 1 provides a quantitative test of the simulation [15 . 

A disadvantage of the present simulation algorithm is that it is not a 
priori clear how many solutions {A{p)i, A{p)*) that Eqs. (13. 9 p and f IS.lOp 
have. Thus we cannot be completely sure whether we have found all the 
solutions or not. Another limitation of the algorithm is that it is applicable 
only to supersymmetric boundary conditions; so it cannot explore physics 
with the finite temperature, for example. 

4. Simulation parameters 

In this paper, we fix the coupling constant in Eq. (11. ip to be 

a\ = 0.3. (4.1) 

This is the same choice of the coupling constant as Ref. 15|, provided that 
our parameter a in Eq. (12. 5p is identified with the lattice spacing in Ref. [15 . 
The size of the momentum grid Nq x A^i is varied from 16 x 16 to 36 x 36. For 
each value of A^o x A^^i, we generated 1280 configurations of the Gaussian ran- 
dom numbers {N (p) , N (p)*) . We then solved Eqs. (13. 9 p and (I3.10p by using 



the Newton-Raphson (NR) method [37|. The solution {A{p),A{p)*) depends 



on the initial guess in the NR method and, as the initial guess, we used 100 
random configurations {A{p),A{p)*) generated with Gaussian random num- 
bers with a unit variant. We judged that the convergence of the NR method 



We would like thank Martin Liischer for bringing our attention to this point. 
^With conventional simulationalgorithms, one needs more elaborate method to com- 
pute the Witten index; see Ref. |3Q] . 



Table 1: Classification of obtained configurations. A is given by Eq. (13.141) and S is defined 
by Eq. ([57 



No = Ni 16 18 20 22 24 26 



+) 

+,+,+) 



6 


1273 


1275 


1271 


1271 


1273 


4 


7 


5 


9 


9 


7 

























































(-,-,+,+,+,+) 




















A 


2 


2 


2 


2 


2 


2 


6[%] 


0.2(2) 


0.1(1) 


-0.2(1) 


0.1(1) 


-0.1(1) 


-0.1(1) 



is achieved when the inaxiinum norm of the residue of Eqs. (13. 9p and (I3.10p 
becomes smaller than 10~^^. Two obtained solutions were regarded identical 
when the the maximum norm of the difference of the solutions is smaller 
than 10^^^. In this way, we obtained configurations tabulated in Tables [1] 
and [21 The total amount of computational time was 34 307.7 core ■ hour on 
the Intel Xeon 2.93GIIz. In Tables [T] and [2], configurations are classified 
according to the number of associated solutions {A{p)i, A{p)*) and the asso- 



ciated sign of the Jacobian [15|. For example, an entry in a row with the 
symbol (— , +, +, +) implies that there was a configuration of {N{p), N{p)*) 
for which we have four different solutions {A{p)i,A{p)*) {i = 1, 2, 3, 4) and 
the sign of the Jacobian is negative for one of them and it is positive for 
other three. Thus, we can read off A in Eq. (I3.14p from this table. We see 
that the expected relation A = 2 for cubic superpotential (11.10 actually holds 
almost within one percent even for the worst case, for A^o x ^i = 36 x 36. 
This fact suggests that even if our root-finding code missed some of solutions 
of Eqs. (13. 9p and (I3.10p . they are precious few. 

5. SUSY WT identities 

WT identities associated with exact symmetries in the formulation pro- 
vide a consistency check for the numerical simulation. In particular, since 
SUSY is an exact symmetry in the present formulation, SUSY WT identities 
must hold even with a finite UV cutoff. This aspect is quite different from 
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Table 2: Classification of obtained configurations (continued). A is given by Eq. (|3.14p 
and S is defined by Eq. 



No = Ni 28 30 32 34 36 



+,+) 1264 1261 1250 1254 1221 

-,+,+,+) 16 17 24 17 26 

+) 1 2 6 31 

+ ,+,+) 4 2 2 

-,+,+,+,+) 1 



(-,-,+,+,+,+) 








1 





A 


2 2.000(1) 


2.002(2) 


1.997(2) 


1.977(4) 


6[%] 


-0.0(1) 0.1(2) 


0.0(2) 


-0.1(3) 


0.0(5) 



conventional lattice formulations in which one hopes that SUSY is restored 
only in the continuum limit. Thus an explicit confirmation of SUSY WT 
identities in the present formulation is interesting in its own right. 

SUSY transformations in the 2D A/" = (2, 2) WZ model are given by Qi, 
Q2, Qi, and Q2- In terms of Fourier modes, they are given by 

QiAip) = -2V2ip,A*{p), QiA*{p) = 0, 

QiF*{p) = 2V2tp-Jiip), Qii'iip) = 0, 

Q^A{p) = V2^,{p), Q,Mp)=0, 

Q,Mp) = V2F{p), Q,F{p) = 0, (5.1) 



for Qi 



Q2^2{P) = -2V2ipJ*{p), Q2A*{p) = 0, 

Q2F*{p) = -2V2ip4i{p), Q2k{p) = 0, 

Q2A{p) = V2^2{p), Q2MP) = 0, 

Q2Mp) = -^F{p), Q2F{p) = 0, (5.2) 



11 



for Qr 



for Q 



i' 



Q.Fip) = -2V2ip-Mp): 
Qii*(p) = v^li(p), 
Qil2(p) = -^F*{p), 

Q2MP) = -2^/2tp,A{p), 
Q^F{p) = 2V2ip,Mp), 
Q2A*{p) = V2i>iip), 
Qiij^ip) = V2F*ip), 



QiMp) = 0, 
QiMp) = 0, 
Qili(p) = 0, 
QiF*{p) = 0, 

QiMp) = 0, 

QiMp) = 0, 

QiMp) = 0) 
Q2F*{p) = 0, 



(5.3) 



(5.4) 



for Qr 



A simple one-point SUSY WT identity is given by 18 



0={S) = {Sb)+(Y,1-2Y,A =(5s)-(iVo + l)(iVi + l), (5.5) 
\ p p / 

where Sb is the bosonic part of the action after integrating over the auxihary 
field, giving by 

Sb^j^Tn*{-p)N{p). (5.6) 



In Eq. (15. 5p . the first equality holds because the action S in Eq. (12. 9p can be 
written as 

S = -\Qyj^ E{ [^WzA*{~p) + F*{-p) + 2W\A)*\Mp) 
° ^ p 

+ \-2ip,A{p) + F{p) + 2W\Ay*\M-p)]. (5.7) 

where Q = — (Qi + Q2)/ ^ is a nilpotent SUSY transformation. Therefore 
(S) = provided SUSY is not spontaneously broken, as the non-zero Witten 
index in the present system shows. One can then exactly evaluate the expec- 
tation value of the parts of the action which are quadratic in the auxiliary 



12 



field and the fermion field; this gives the second equality in Eq. fIS.Sp . Since 
J2p 1 = (^0 + l)(^i + 1)) we have the last equahty in Eq. fl5.5p . That is, we 
have 

This relation provides a quantitative test of the simulation 18|, llSJjj We 
show 6 in Tables [Hand [2] for each value of A^o ^ ^i and see that relation (I5.8p 
holds within 0.5%. 

As for a two-point SUSY WT identity, from the relation {Qi{A{p)ijJi{—p))) 
0, we consider 



Po (^A{p)A*{-p)^ = - Im {jijji{p)ipi{-p)J 



(5.9) 



In Fig. [H we plotted both sides of this relation for our coarsest grid ( A^o x A''i = 
16 X 16) as a function of apo; both sides coincide within the statistical error. 
For more quantitative comparison, we plotted in Fig. |2] the ratio of the 
left-hand side and the right-hand side again as a function of apo- (The error 
in the ratio was estimated by a simple propagation of error rule.) We plotted 
the cases of A'o x A?"! = 16 x 16 and A'o x A?"! = 36 x 36. For both cases of the 
momentum grid, the identity holds within ~ 5%, a good indication of SUSY. 
Similarly, from {Q2{F*{p)i'i{—p))) = 0, we havqj 

F{p)F*{-p) 



-- -piRe {i'i{p)i'i{-p)j + polm.{i)i{p)i)i{-p)) . (5.11) 
The results for this identity are depicted in Figs. |3] and Ho 



^Lattice formulations adopted in Refs. llSll and [ISl and more generally lattice formu- 
lations based on the Nicolai map 23|, [ITI, llSl, |l6|, |2J, [19, 20, 21j] possess exact invariance 
under the above nilpotent Q transformation. See also Ref. [38J and references therein. 

^In the Wess-Zumino model, the integration over the auxiliary field is defined by an 
analytic continuation, and according to this prescription, the correlation function is ex- 
pressed as 



F{p)F*{-p)) 



Nip) + 2ip,A{p) -N*{-p)-2ipiA*{-p) 



LnL 



O-i^l, 



(5.10) 



in terms of a correlation function of the Gaussian random fields and the scalar fields. The 
last term gives a negative contribution, and despite its appearance, the left-hand side can 
be negative as our plot shows. 

^Since (F(0)i?*(0)) = -(Q2(V'i(0)i?*(0)))/%/2, {F{p)F*{~p)) evaluated at p = must 
vanish if SUSY is not spontaneously broken. Our plot in Fig. |3] is consistent with this 
expectation. 
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Figure 1: The left-hand side ("bosonic") and the right-hand side ( "fermionic" ) of Eq. (|5.9I 
as a function of apo along the line api — 0; Nq x Ni = IQ x 16. 



Having obtained these encouraging results concerning SUSY in our nu- 
merical simulation, let us proceed to study physical questions. 

6. Conformal dimension of a chiral primary field 

As noted in the Introduction, in the IR limit, the scalar field A in the WZ 
model with cubic potential f ll.ip is expected to behave as a chiral primary 
field with the conformal dimensions {h, h) = (1/6, 1/6). For such a field, the 
two-point function will behave as 



{A{x)A*{0)) oc 



^2h^2h ' 



for |x| large. 



(6.1) 



Then assuming that A is spinless, h = h, the anomalous dimension h+h = 2h 
could be extracted from the susceptibility of the scalar field x^, defined by [15 



X<t> 



a JlqLi 



d^x{A{x)A*{0)) 



o^LqLi 



m 



(6.2) 
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Figure 2: The ratio of the left-hand side and the right-hand side of Eq. (|5.9p as a function 
of apo along the line api — (the origin p = is excluded from the plot). The cases 
of A^o X TVi = 16 X 16 and A'^o x TVi = 36 x 36 are plotted. 



as 



X0 oc 



(Px ■ 



1 



LqLi 



(X 



2\2h 



OC (LqLi 



\l-h-h 



for Lii_ large. 



(6.3) 



Thus \n{x^) would be a linear function of In(LnLi) for large L^ and 1 — h — h 
is obtained by the slope of the linear function jlSj. In Fig. El we plot ln(X(^) 
as a function of ln(a~^LoLi) (recall that in our present simulation the lattice 
spacing is fixed to aX = 0.3). From the plot, we see that a fit by a linear 
function would be good for large sizes as q'^LqLi > 24 x 24. Therefore, we 
applied a linear x^ fit by using data from A^q x ^i = 24 x 24 to A^o x ^i = 
36 X 36. To estimate a systematic error associated with this choice of fitting 
range, we also carried out a linear fit using data from A'"o x A^^i = 26 x 26 
to A^o X A^i = 36 X 36. See Tabled 

In Eq. (16. 2p . the integral over x is performed for all x including the co- 
incidence point X = 0. This might be physically unnatural because the 
coincidence point x = would suffer from ambiguity associated with the 
UV regularization. On the other hand, true long-distance physics should be 
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Figure 3: The left-hand side ("bosonic") and the right-hand side ( "fermionic" ) of Eq. (|5.1ip 
along the line api = 0; No x Ni = 16 x 16. 



independent of such a UV ambiguity. To avoid such UV ambiguity and to 
estimate how much our determination depends on the UV prescription, we 
could define the susceptibihty with the UV part subtracted [15|: 



d^x{A{x)A*{0)) 



^ J LQL\—non\o? 



u'^LqLi 



i(0) 



a?{LQL 



1 v^ 2 . fpoano\ 2 . fpiani\ 
^ ,„ > — sm — sm 



A{p) 



(6.4) 



where the integral has been defined by extracting a region a^riQ x ni con- 
taining the coincidence point x = 0; in this expression it is understood that 
(2/p^) sin(p^an^/2) = an^ for p^ = 0. Figure [6] is the result with the UV 
subtraction with ng x ni = 3 x 3 (this is a choice identical to Ref. 15|). 
Our linear fits are summarized in Table |3J By adopting the average of the 
second and the third rows in Table [3] as the central value and estimating a 
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Figure 4: The ratio of the left-hand side and the right-hand side of Eq. (|5.1ip as a function 
of apo along the line api — (the origin p = is excluded from the plot). The cases 
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systematic error associated with the fitting range and the UV ambiguity by 
variations in Table |31 we quote 



1 -/i-/i = 0.616(25)(13). 



(6.5) 



This value is somewhat smaller than the value obtained in Ref. [15|, but 
still consistent with the expected exact value 1 — h — h = 2/3 = 0.666 . . . 
within 1.3(7. 



7. Central charge from a supercurrent correlator 

The central charge c, defined by the Virasoro anomaly 

1 c/2 



(T,,(a;)T,,(0)) 



^7r2 z^ 



(7.1) 



where Tzz is the holomorphic part of the energy-momentum tensor, is a fun- 
damental quantity that characterizes a conformal field theory (CFT). It is 
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Figure 5: ln(x0) as a function of ln(a~^Loii). No UV subtraction is made. The broken 
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thus of great interest whether we can measure this quantity c in the IR re- 
gion from the present numerical simulation and obtain further support for 
the conjectured correspondence between the WZ model and SCFT. A salient 
feature of the present formulation is its nice symmetry properties, including 
translational invariance. Thus, we may define the energy-momentum ten- 
sor as a Noether current associated with the translational invariance in a 
regularized framework. 

In an A/" = 2 SCFT, the central charge c appears also in the correlation 
function between the holomorphic U{1) currents: 



[Ux)JM) 



1 c/3 

^^2 ~;2 



(7.2) 



In the LG description, the U{1) symmetry in the M = 2 SCFT is identified 
with the U{1)r symmetry of the WZ model with a (quasi-) homogeneous 
superpotential (see below) and this U{1)r is also exactly preserved in the 
present formulation. We can thus define a conserved U{1)r current whose 
appropriate component is identified with J^. 



Table 3: Linear ^ fit used for a determination of tlie anomalous dimension in Eq. (|6.5I 
no X n\ denotes the UV subtraction region defined in Eq. (|6.4p . 



rio X nx 


Fitting range of A^o x A^i 


xVd.o.f. 


\-h-h 


0x0 


from 24 x 24 to 36 x 36 


0.904 


0.603(19) 


0x0 


from 26 x 26 to 36 x 36 


1.088 


0.609(25) 


3x3 


from 24 x 24 to 36 x 36 


0.910 


0.624(20) 


3x3 


from 26 x 26 to 36 x 36 


1.108 


0.629(26) 



In preliminary numerical study of the above two-point functions of bosonic 
operators {Tzz for Eq. (17. ip and Jz for Eq. (17. 2p ). however, we found that 
signals are generally very noisy and reliable fits for the central charge would 
be rather difficult]^ So, in what follows, we adopt a different approach that 
employs a two-point function of the M = 2 superconformal currents: 

(G+(a:)G;(0)> = -i^^. (7.3) 

Since the operators G^ are fermionic, there is no "disconnected diagram" 
which could contribute to Eq. (17. 3p and we expect a clear signal from numer- 
ical simulation. In fact the following results seem to be in accord with this 
naive expectation. 

Now, since our present formulation possesses exact SUSY, we can define 
a conserved supercurrent. To obtain this, we consider the localized SUSY 
transformation (for a generic field (p) 



^^(P) = L^X 



Y.~^''^q)Q.(p{p - q) -J2tiq)QaHP ' q) 



, (7.4) 



where ^"(g) and ^"(g) are non-constant Grassmann parameters and Qa 
and Qa are defined by Eqs. (I5.ip -f l57i|) . Then from the variation of the 



^We stress that this preliminary study was quite incomplete. In particular, we have not 
seriously considered the necessity of subtracting some disconnected parts from correlation 
functions. Clearly additional study is needed to conclude something about the use of those 
correlation functions in computation of the central charge. 
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action, we define supercurrents 5'^(p) and 5'^(p) as 

(7.5) 
By construction, above supercurrents satisfy SUSY WT identities, such as 



Pm (S^iP)'PliQl) ■ ■ ■ ^niqn)) = 7;Y1 (^l(^l) • • • Q2Vi(9i + P) ■ ■ ■ Vn{qn)) ■ 



i=l 



(7.6) 

This identity would imply that (according to the standard argument) the 
current S'^ is a correctly normalized operator. 
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In Eq. 07. Sp . the explicit form of the supercurrents are given by 

Stip) = tV E ^'(P - ^)^^(P - ^)^2(g), (7.7) 

stip) = TT^J2 -^w'{A){p - q)Mq), (7.8) 

S;{p) = --^ J2 ^'^P - ^)-^*(P - ^)^2(g), (7.9) 

~sj{p) = yt;Y. -^^'i^Tip - ?)^i(9), (7.10) 

q 



and 



§^ip) = --^ V y2^(p - q),A*{p - q)^P,{q), (7.13) 

^;(P) = -;^ E 7f ^'(^)*(^ - ^)^2(g). (7.14) 



A Noether current is, however, always ambiguous as one can change the 
definition as Sf{p) — )■ Sf{p) + p^X^ip) and Sf{p) — )• Sf{p) — p^X^ip) by 
using certain combinations X'^{p) without affecting the conservation law, 
PzSf{p) + PzStip) =0. To fix this ambiguity, we required in Eqs. (17. 7p 
to (EIID that 

St = S^ = Bt = k = 0, (7.15) 

when W = 0. The WZ model with W = is a massless free theory that itself 
is an A/" = (2, 2) SCFT. For this system it is natural to adopt a supercurrent 
that obeys the gamma-traceless condition. 



E^mS =E^M?- =0' (^-16) 



M / /x \ M 
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because this condition is a super-partner of the traceless condition T^ 



[39|] that is usually assumed in CFT. Eq. ( 17.1 5p is nothing but Eq. ( I7.16P 
in components. 

Thus we had a natural definition of the supercurrents. We then postulate 
a correspondence between the components of above supercurrents and the 
holomorphic part of the superconformal currents in the IR limit. That is, 

Si ip) ^ Gt ip) , S; ip) ^ G; ip) . (7. 17) 

Our reasoning for this correspondence is as follows: 

(i) The conservation law of the supercurrents in the coordinate space yields 
dzSf{x) + dzSf{x) = 0. In the IR limit, the derivative of the superpo- 
tential iy(74) is expected to become irrelevant (see, for example. Sec- 
tion 14.4 of Ref. O) so Sf{x), from explicit forms ([7SD and fmojl . 
could be neglected in the IR limit. Then the conservation law implies 
that Sf{x) are holomorphic functions, a correct property of G'^(x). 

(ii) The WZ model with a homogeneous superpotential W = A$"/n pos- 
sesses the U{1)r symmetry 

i,^ ^ e(^-2/")*'^^<i, F ^ e-(2-2/«)^^F, (7.18) 

which is, in the IR, identified with a f/(l) symmetry generated by a sum 
of two f/(l) charges (g, q) in the left- and right-moving M = 2 super- 
conformal algebras (see, for example. Section 19.4 of Ref. [ll|). Thus, 
in Eq. (I7.18p . the U(1)r charge of the scalar field A is assigned 2/rz, 
because A is identified with a chiral primary field in the A^-i model 
with the U{1) charges (g, g) = (1/n, 1/n). Then, assuming that S~i 
in Eq. (17. 7p possesses U{1) charges (g, 0) in the IR, we see from the 
t/(l) charge assignment in Eq. (I7.18P that q = +1, a correct U{1) 
charge of G^. Similarly, for S~ in Eq. (17. 9p . we have q = —1, a correct 
U{1) charge of G~. 

(iii) The overall normalization of the supercurrents has been fixed such that 
Eq. (17. 3p is reproduced under correspondence (I7.17P for the massless 
free theory [W = 0). This theory itself is an A/" = (2, 2) SOFT whose 
left-moving sector possesses the central charge c = 3 (one complex 
scalar and two Majorana-Weyl fermions). 
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Thus, under identification (I7.17p . our procedure is as follows: we numerically 
compute the two-point function of the supercurrents in the momentum space 

{snP)S;{-p)). (7.19) 

Then, in view of correspondence (17.17p . we compare this function in the IR 
regioio |p| ^ -^ with the two-point function of the superconformal currents 
in the momentum space 



LoLi f d^xe-'P^{GUx)G;iO)) 



- ^oL.^|;f;i..(b|5) ^^ LoL,-f ^, (7.20) 

487rc/p|d)^ 24:71 pz 

where we have used Eq. (17. 3p that is, strictly speaking, an expression valid 
only on M^. In deriving the last expression, we regularized the singularity in 
the integrand at x = by setting 1/z^ = {z)^/{x'^)^ — )■ {z)^/{x'^ + 6"^)^. Note 
that the last low- momentum behavior (i.e., the IR physics) in Eq. fl7.20p is 
independent of the regularization parameter S. 

For the above computation, we used only data with the finest grid Nq x 
A^"! = 36x36. In Fig. [71 we depicted the real part of correlation function (17.190 
as a function of apo along the line api = vr/lS ~ 0.1745. The broken line is 
the real part of function (17.200 with c = 1; it agrees well with the data for 
\o-Po\ < vr/18 ~ 0.1745 < 0.3, a good indication for c ~ 1 in the IR region. 
Fig. [8] is the same as Fig. [3, but for the imaginary part. The broken line is 
the imaginary part of function (17.200 with c = 1; this time it agrees well with 
the data for |apo| < 27r/18 ~ 0.349 ~ 0.3. 

Our fit for c proceeds as follows: We first consider the real part of the ratio 
of correlation function (17.190 to the function LoLi{i/24:7i)pl/p2 appearing 
in Eq. (I7.20p . For illustration, the real part of this ratio along several constant 
api lines are depicted in Figs. [H [10] and [Til We see that actually the ratio 
is close to 1 around the origin in the momentum space \p\ ~ 0. Then, to 
use data points with a fixed energy scale, we define a fitting region R{bi, 62) 
on the momentum grid as illustrated in Fig. [121 Then we carry out a x^ fit 
by a constant for data points in the fitting region R{bi,b2)- Since the IR 



^''Notc that the coupling constant A in Eq. (|l.ip is a unique physical dimensionful 
parameter in the present system on M^. 
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region is characterized by the condition \ap\ ^ aA = 0.3 and the lowest 
non-zero momentum in the present iVo x iVi = 36 x 36 momentum grid is 
ap^ = 7r/18 ~ 0.1745, we regard the fitting region i?(l, 1), or at most -R(l, 2), 
as the IR region; see Table HI From the numbers in the table, we estimate 



c= 1.09(14)(31), 



(7.21) 



as the central charge in the IR region. Our estimate reproduces the conjec- 
tured value c = 1 very well. 



Table 4: Results of a x^ fit used for a determination of the central charge c in (|7.21l) . 



Fitting region Number of data points x^/d.o.f. c 



i?(l,l) 
i?(l,2) 



8 1.93 
24 3.50 



1.09(14) 
1.40(8) 
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Figure 8: The imaginary part of the correlation function (|7.19p as a function of apQ along 
the line api = 7r/18 ^ 0.1745; A'o x A^i = 36 x 36. The broken line is the imaginary part 
of function ([7:201) with c = 1. 



Although strictly speaking expression f l7.20p is meaningful only for very 
small momenta and the LG description is expected to be valid only for the 
low-energy region \ap\ ^ aA = 0.3, we found that it is nevertheless interesting 
to repeat the above procedure for "intermediate" energy \ap\ > aA = 0.3 or 
even "high" energy \ap\ ~ tt regions. Figures [T^ and [T^ is a result of the 
fitting in an "intermediate" energy region \ap\ ~ 0.7. We observe a rather 
good fit with an "effective c", c ~ 2. We may repeat such a fit by changing 
the parameter h in R{h, b) which roughly corresponds to the energy scale. 
The result of such fits as a function of b are shown is Fig. [151 It is interesting 
that, going from the UV side to the IR side, the plot goes from c ~ 3, 
the central charge corresponding to the massless free theory, and decreases 
to c ~ 1, the central charge of the A2 model. Although one should not take 
this plot so seriously because x^ becomes quite large for large b (for example, 
X^/d.o.f. ~ 19 for 6 = 8), it is still interesting to note an analogue of the 
Zamolodchikov C function [40(] in this plot; the C function monotonically 
decreases along the renormalization group flow from UV to IR and coincides 
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Figure 9: The real part of the ratio of the correlation function (|7.19p and the func- 
tion LoLi{i/24Tr)pl/pz as a function of apo along the line api = 0; Nq x Ni = 36 x 36. 
The origin p = is excluded from the plot because pl/ps is singular at the origin. 



with the central charge at a fixed point. It is an interesting problem to clarify 
the underlying reason for this similarity of Fig. [15] with the C function. 



8. Conclusion 

In this paper, we carried out a non-perturbative numerical study of low- 
energy physics of the 2D WZ model with the massless cubic superpotential. 
We obtained the critical exponent 1 — h — h and the central charge c which 
are consistent with the conjectured emergence of a non-trivial Af = (2, 2) 
SCFT, thus provided further support for the LG description. Our results 
indicate that our supersymmetric non-perturbative formulation of the WZ 
model is working, although there has been an issue concerning the locality in 
this formulation. For further improvement of these results, better observables 
which yield less systematic errors should be investigated. Also, we want to 
generalize the present analysis to higher critical LG models. 

We would like to thank Michael G. Endres, Kazuo Fujikawa, Masafumi 
Fukuma, Daisuke Kadoh, Hikaru Kawai, Akitsugu Miwa, and Tsuneo Ue- 
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